An analytically tractable model of bad metals 
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We discuss a model Kondo-type Hamiltonian representing an analytically tractable version of the 
model used by Yin et.ai, Phys. Rev. B86, 2399 (2012) to explain the non- Fermi liquid behavior of 
iron chalcogenides and ruthenates in an intermediate energy range. We consider a regime where a 
complete screening of the local degrees of freedom proceeds in two stages described by two charac- 
teristic energy scales Tf[ b >> Eg. The first scale marks a screening of the orbital degrees of freedom 
and the second one marks a crossover to the regime with coherent propagation of quasiparticles. 
We present analytical results for the specific heat and magnetic susceptibility at T « T^[ h . 



I. INTRODUCTION 



\^ ' There is a significant number of metallic systems dubbed "bad metals" where the energies at which quasiparticles 
emerge as coherent objects are much lower than the characteristic scale of the interactions. In [H,[II it was suggested 
that in the compounds based on d and /-elements, the Hund's interaction plays major role in delaying the onset of 
coherence. Those authors considered models where electrons carry both orbital and spin indices using the Dynamical 
D ■ Mean Field Theory (DMFT) and found that the coherence scale Eq was indeed low in comparison with the bandwidth 
W or the Hund's rule coupling Jh- In the intermediate range Eq < T,u <min(VF, Jh) the electronic self energy was 
found to have a non- Fermi liquid form £ ~ u b . In [l[ the exponent was found to be universal b = 1/2, in [2| which 
used a more realistic model b was non-universal. The latter result was found to fit the observed behavior of the 
mid-infrared optical conductivity in the iron chalcogenides and ruthenates. 

The starting point for 2] is the Kondo lattice model where the localized <i-electrons (6 electrons per site) give rise 
I , to spin and orbital moments: 

g ■ H K = <Wt (fc)^v(fc) + (1) 



where N is the number of sites. The spin degrees of freedom are described by spin S = 2 operators acting on the spin 
indices and the orbital sector is described by the Hubbard operators Xji (j, 1 =1,...M). The symmetry of Hamiltonian 
(HJ) is SU(M)xSU(2)xU(l). In iron pnictides and chalcogenides, a Fe-ion is surrounded by a tetrahedron of pnictogen 
\ or chalcogen, and the resulting crystal field is weak in comparison to the Fe-pnictogen hybridization. As a result all d 
■ orbitals have degeneracy, M = 5. In the ruthenates where the coordination of Ru-ion is octohedral the crystal field is 
strong, yielding M = 3 because only the ti g orbital actively participates in the interaction. A detailed description of 
£f*^ . the model can be found in the Supplementary material to [2J . Furthermore, a similar model with different parameters 
t-H ' was considered in [3] . 

As we have mentioned above, model (|T|) was treated by DMFT 0. In that approach local moments from different 
• i-H ' sites have been considered as uncorrelated. As a result the Kondo lattice was treated as a single impurity problem 
with a self-consistently renormalized density of states (DOS) of the band electrons. It turned out, however, that in the 
given case the DOS at the chemical potential remains non-singular and hence the self-consistency had no qualitative 
. . .' effect. Therefore the range of energies where the magnetic interactions between sites are still small can be treated as 
a single impurity problem. 

Additionally, the DMFT treatment established the existence of the intermediate regime marked by nontrivial power 
laws in the electron self energy. In all likelyhood this regime emerges as a crossover between the quantum critical 
point (QCP) of the purely orbital Kondo model (Ji = 0) and the Fermi liquid strong coupling regime of the full 
model (J2 7^ 0). Indeed, at J2 — the single impurity Kondo model decouples into two independent Kondo models 
describing scattering in the orbital and the spin sectors. This decoupling occurs due to the fact that electronic 
densities ("currents") coupled to the orbital Xji and spin operators S a commute with each other. Both spin and 
orbital Kondo models are overscreened and at strong Hund's coupling, the sign of the exchange interaction in the 
spin channel is ferromagnetic 0. Therefore, at J2 = the orbital Kondo model scales to the QCP characterized 
by nontrivial exponents and the spin one scales to weak coupling. However, as soon as the coupling J2 is switched 
on, the renormalization group trajectories start to deviate from the critical point. Eventually the system becomes a 
Fermi liquid. The corresponding energy scale is, of course, determined by J2. 
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Although the meaning of the numerical results by [2| is sufficiently transparent, it would be advantageous to have 
a simple model of bad metals where a qualitatively similar picture can be obtained analytically. Therefore, this is the 
primary focus of the present paper, which is organized as follows: The model is described in Section II; in the same 
section we present its solution and describe the low temperature thermodynamics. In Section III we briefly discuss 
the general case (that is the single impurity version of model (JTJ). Lastly, we provide an overview of our results in 
the conclusion. 



II. A SOLVABLE SINGLE IMPURITY MODEL 



Below we consider a single impurity version of model ([T]) where an analytical treatment is possible. This version 
has four species of fermions with symmetry SU(2)xSU(2)xU(l) interacting with pseudospin T and spin S, both of 
magnitude 1/2. The interaction has a form: 

V = gi^ + (I ® r a )^f a + 2g 2 4> + {a a ® T b )ipS a f b + g 3 ^ + {a a <g> I)^S a (2) 

The Pauli matrices a a and r a act in the spin and the orbital sectors respectively, and obey the identites, 

a a S a -1/2 = P spin , r a f a + 1/2 = P orb , (3) 

where P sp i n , P rb are permutation operators acting in the corresponding sectors. Therefore, when gx = g% = .93 
the interactions become P sp inPorb — Psu(4)> which yields the integrable SU(4) Coqblin-Schrieffer modelQ- At low 
temperatures the latter model displays Fermi liquid behavior characterized by the complete screening of the impurity. 
When the couplings are not equal to each other the model is not integrable (except for the case g% = 0), but in all 
likelyhood eventually reaches the regime of full screening. 

We will be interested in the case <?i >> (72,3 when the orbital sector reaches the critical point first. As is well 
known, for the single impurity problem, the dimensionality of the bulk is irrelevant as long as the density of states 
(DOS) at the chemical potential is constant. It allows one to treat the bulk in the single impurity Kondo problem 
as a one-dimensional theory of chiral fermions. Such a replacement carries many advantages enabling one to apply 
various non-perturbative techniques available for one-dimensional theories, such as Bethe ansatz and non-Abelian 
bosonization described in [5], [6]. In the latter procedure one replaces the bulk Lagrangian of ID fermions with 
U(l)xSU(2)xSU(2) symmetry by an equivalent representation consisting of a U(l) bosonic theory describing charge 
fluctuations and two SU2(2) Wess-Zumino-Novikov-Witten (WZNW) Lagrangians describing the spin and the orbital 
sectors. The peculiarity of the N — M — 2 case is that the WZNW models are equivalent to models of Majorana 
fermions. The net result is: 

J2^U k )( d r-^a(k)= (4) 

k 

f°° r 1 1 

/ dx d x 4>(id T + d x )(j) + -Xa{d T - id x ) Xa + ^£,a{d T ~ id x )£ a 

where £ a ,Xa (a = 1,2,3) are Majorana fermions transforming according to the adjoint representation of the SU(2) 
group. The scalar bosonic field 4> describes the charge sector which is decoupled from the impurity. Henceforth, we 
set the Fermi velocity vp — 1 and thus p(ep) = l/2n. The total central charges of the left- and right-hand side of (|4j) 
are, naturally, equal: 4= 1 + 3/2 + 3/2. It follows that the electron current operators can be expressed as 

tA+tV -> \e abc XbXci ^+a a i> ^e abc te c (5) 

^ + rVV -> i X a t b (6) 
As was shown in at the 2-channel Kondo model QCP the local pseudospin T a renormalizes into [§| 

f a ^T° K rb - 1/2 e Xa (0) (7) 

where e is a local zero energy Majorana mode residing on the impurity site x = and T^f b ~ Wg\ exp(— n/gi) is the 
orbital Kondo temperature. 

As a result we obtain the following effective Lagrangian describing the behavior at energies below Tj[ (the U(l) 
part of the bulk Lagrangian is omitted): 

L = gUT£ b ) i/2 ie(o)s b }e+i^a b cS a em c (o) + ^ I 



\ f dx Xa {d T ~id x ) Xa + l -ed T e + T° K T b 1/2 e X M X2 {Q) X M . IM 
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where g^g^ are renormalized values of the corresponding coupling constants. The term in the square brackets is the 
critical point Lagrangian of the orbital 2-channel Kondo model. The last operator in the square brackets is irrelevant, 
but it is included since it determines the impurity contribution to the specific heat @: 

C v ~ {T/T° K rb ) ]n(T£ b /T). (9) 

The most important part of (J5J is the first term, which follows from the gi term in Eq.@ where we replaced the 
fermionic bilinear by (j6j), the orbital operator T a by (J7J and applied the fusion rule: 

92 



2g 2 ^ + (a a ® r b )^S a f b -> g* 2 (T%~ b )- 1/2 {x b C}(t + V )[S a X b }(t)e(t) a [?§ a ]e (10) 

port 
K 



T° rb 2TTT] 



Since the interaction becomes retarded in the process of renormalization, we have assumed that the operators should 
be time split by the amount 77 ~ 1/T^' b . As a result, at the QCP the operator in question becomes a relevant one 

To solve model © we introduce a transformation 

f a = 2ieS a , (11) 
where f a are Majorana fermions, satisfying the Clifford algebra 

[f a J b ] + = 25 ab . (12) 
This transformation respects the Casimir operator S 2 = 3/4. Since e — 1, the inverse transformation is 

S a = ^abcf b f c , e = /i/a/ 8 . (13) 

The latter expression reproduces the commutation relations of spin operators. It is essential that / a 's replace com- 
pletely e, S a and therefore field e is no longer in use. Then the effective Lagrangian ((HJ becomes 

L = i9* 2 ^/rf b [e(0)f b ] + \j dxUdr - id x )U + \f a d T f a (14) 

+ [\J dx Xa (d T -id x )xa + f [/V(0)][/T(0)] +T° rf ~ 1/2 /i/ 2 /3Xi(0)x 2 (0)x3(0)]. (15) 

The impurity thermodynamics is determined by the Green's functions of the /-operators. In the zeroeth order in g% 

1/0 

and Tf[ b , one obtains these functions by diagonalizing the quadratic part of the Lagrangian given by ([M)) . The 
most convenient way to proceed is to write down the Lagrangian for the Fourier components of / and £(0): 

L ° = E [^H"' °) G o V> * = 0)£>, 0) + iv/Eb72£ a (-")/° M + y / a (-")/»] (16) 

where 

G (w,a; = 0) = / ^^-7: = -^sign(w) 
J 2tt iui — k 2 

is the Green's function of the bulk fermion at x — and Eq = \\g2^T^ h ' . The net result is 



( «/M/(-w)» ((/Mg(-C,0))) \_ 1 /isigl^n -i^72\ r , 

^((^,0)/^))) ((^,0)C(-c,0)»; |w n | + £b ViV^72 i^n/2 7' 1 J 

From the form of (TlTl) is clear that E'o is the crossover energy scale to the Fermi liquid regime. 

The nonquadratic terms in (|15[) provide subleading corrections to the scaling and hence are irrelevant in the infrared. 
The corresponding corrections can be calculated by the perturbation theory using Green's functions (|17j) . For instance, 
the ferromagnetic exchange interaction 53 reduces the value of Eq so that the new scale E e ff is given by: 

E%) = El 12 - 2 5 3V / 2iT^(/(u,)£(-u,,0)> ) 



_gl/2 _ / 18 \ 

e// ~ 1 - 5a A HTfrt/maxfT, E ef f )]' V > 
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FIG. 1: Temperature dependence of the magnetic susceptibility ([20 



The onset of Fermi liquid is delayed till the new scale E* determined by T = limit of equation (|18l) : 

E*[l - ^ ln{T£ rb /E*)} 2 = E 



(19) 



which can be significantly smaller than Eq if the ferromagnetic exchange is strong. Otherwise the weak temperature 
dependence of E e fj will give logarithmic corrections to the scaling considered below. In the first approximation the 
thermodynamics is still determined by (|20I22|) with Eq replaced by E e ff(T) determined by fp~8]) . The last term in 
(I15p containing six fermionic operators is even less singular at T « T^ b . 

In the first approximation one can neglect the temperature dependence of E e ff given by (|18[) and replace it by Eq. 
The temperature dependence of the magnetic susceptibility (displayed on Fig. 1) is given by the formula [Til ] 



« S Z S Z »= x (T,w = 0) = - 



UjEn 



nE 



1 - (2ttT/E ) 



irj {^+EtY 
- 1 ^'(£ / 27 rT+l/2 



tanh(w/2T)dw 



which at T > O.SEq can be approximated with high accuracy as 

X(T) ' 



The spin sector contribution to the specific heat is 
3T r dxx 2 

O — 



4(T + 0.59^0) 



3T 



4irE J [x 2 {T/E ) 2 + 1] cosh'(a;/2) 47r£ 



1 



Eo 
2ttT' 



iP'(e /2ttT+1/2 



(20) 



(21) 



(22) 



We see that at T << Eq the susceptibility is constant and the specific heat is linear in T. Hence at energies smaller 
than Eq the impurity spin is fully screened. 



A. Spin-orbit coupling 



Our model ^ can be easily augmented by inclusion of the spin-orbit coupling: 

H so = ATS. 
At the 2-channel QCP we replace T as in ([7]) to get 

Hso -+ i^^e X a S a 

rporb 



(23) 



(24) 
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FIG. 2: Temperature dependence of the spin sector contribution to the specific heat ([22 



This term combines with (|10p such that now the impurity spin interacts with a linear combination of the orbital and 
spin Majoranas: 



cosax a + Binaf, tana = g* 2 T% h /k, 

with the coupling constant 

V = {A 2 /T£ b + {g* 2 ) 2 T1[ b ) 1 ' 2 . 
Thus, we see that the spin-orbit coupling does not bring any qualitative changes. 



(25) 
(26) 



III. A MODEL WITH GENERIC SYMMETRY 



As mentioned earlier, a single impurity version of model ([T]) is not solvable. However, it is still possible and 
worthwhile to make some qualitative statements about it. Consider model (JTJ with U(l)xSU(N)xSU(M) symmetry. 
An analog of decomposition (j4]) in this case is (see, for example, 



M N 
j=l a =l k 

27T 



(27) 



dxd x <j)(id T + d x 



N + M 



dx 



J a J a 



papa 



where J a (a = 1, ...N 2 - 1) and F a (a = 1, ...M 2 - 1) are currents from SU(N) M and SU(M)jv Kac-Moody algebras 
defined as 



M N 
J'=l 



(28) 



(7 = 1 



with cr a , r a being generators of the corresponding groups. For general N =/= M the quadratic forms of currents cannot 
be expressed as models of free fermions as it happens for TV = M, where they can be written in terms of Majorana 
fermions transforming in the adjoint representation of the SU(N) group. Neither does it work for the cross term in 
@. Instead we have 



(29) 



where $, A are Wess-Zumino primary fields from the adjoint representations of the SU(M) and SU(N) groups respec- 
tively. Their conformal dimensions are 



M 
N + M' 



N 
N + M' 



(30) 
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The QCP corresponds to the overscreening of the SU(N) sector. The corresponding Kondo scale is T£ rb — 

WgY^ N exp(— 2ir/Ngi). At the QCP one operator $ creates an average with the impurity spin leaving behind a 
zero mode, fashioned along the lines of (|7|10[) . As a result, the operator perturbing the critical SUjv(M) WZNW 
model is A b coupled to the SU(M) impurity spin. Its scaling dimension is N/(N + M) < 1 and hence it is relevant. 
Thus, we can be confident that the model scales to strong coupling, which is presumably a Fermi liquid, but the 
details of this are not known. Consequently, the corresponding energy scale marking the crossover to Fermi liquid is 

Eo ~ gl ,{1 ~ hK) T% b (31) 



IV. CONCLUSION AND ACKNOWLEDGEMENTS 



We have presented an analytically tractable single impurity Kondo model where the impurity carries both spin 
and pseudospin S = T =1/2 and the bulk electrons carry both spin and orbital indices. Our original motivation 
for considering this model was its relevance to the problem of "bad" metals, as formulated in 0,0], 0- As we have 
already mentioned, according to DMFT, the single impurity problem describes both thermodynamics and the electron 
self energy within a range where the intersite correlations are relatively weak. Although the experimentally relevant 
models have different values of S, we have chosen S = 1/2 to obtain analytic results. 

Furthermore, we have considered the most interesting regime when the orbital moment is screened first and the 
coupling between spin and orbital channels is weak. Therefore, the screening of the orbital moment leads to a non- 
Fermi-liquid quantum critical point (QCP) which is destabilized by the interaction term mixing spin and orbital 
channels. If the corresponding coupling constant 172 is small, there is a wide crossover regime between the unstable 
two-channel Kondo model QCP at high energies and coherent Fermi liquid at low energies. The ratio of these two 
energy scales is ~ g%. 

It is quite likely that the nonuniversal exponents in the self energy Q discussed in the Introduction are crossover 
effects. To illustrate this point one needs to calculate the self energy; we plan to do it in a subsequent paper. 

AMT is grateful to G. Kotliar, S. Lukyanov, R. M. Konik and P. Coleman for inspirational discussions and to A. 
James for help with the numerical calculations. The work was supported by the US DOE under contract number 
DE-AC02-98 CH 10886. 
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